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Abstract 

The  Pallet  Loading  Problem  (PLP)  maximizes  the  number  of  identical  rectangular  boxes  placed  within  a  rectangular 
pallet.  Boxes  may  be  rotated  90°  so  long  as  they  are  packed  with  edges  parallel  to  the  pallet’s  edges,  i.e.,  in  an  orthogonal 
packing.  This  paper  defines  the  Minimum  Size  Instance  (MSI)  of  an  equivalence  class  of  PLP,  and  shows  that  every  class 
has  one  and  only  one  MSI.  We  develop  bounds  on  the  dimensions  of  box  and  pallet  for  the  MSI  of  any  class.  Applying  our 
new  bounds  on  MSI  dimensions,  we  present  an  algorithm  for  MSI  generation  and  use  it  to  enumerate  all  3,080,730  equiv¬ 
alence  classes  with  an  area  ratio  (pallet  area  divided  by  box  area)  smaller  than  101  boxes.  Previous  work  only  provides 
bounds  on  the  ratio  of  box  dimensions  and  only  considers  a  subset  of  all  classes  presented  here. 

©  2006  Elsevier  B.V.  All  rights  reserved. 
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1.  Introduction 

We  identify  each  instance  of  a  Pallet  Loading  Problem  (PLP)  by  a  quadruple  ( X,Y,a,b ).  We  have  a 
rectangular  pallet  with  length  X  and  width  Y  (X  ^  Y),  and  a  rectangular  box  with  length  a  and  width  b 
(a  ^  b ).  Boxes  may  be  rotated  90°  so  long  as  they  are  placed  with  edges  parallel  to  the  pallet’s  edges,  i.e., 
the  packing  must  be  orthogonal.  We  can  assume,  without  loss  of  generality,  that  X,  Y,  a,  b  are  positive  integers 
(e.g.,  Bischoff  and  Dowsland,  1982).  We  also  assume  that  at  least  one  box  can  be  packed  in  the  pallet:  X  ^  a 
and  Y  ^  b. 

We  encounter  PLP  when  trying  to  maximize  the  number  of  identical  boxes  with  dimensions  a  and  b,  placed 
on  a  pallet  with  dimensions  X  and  F  where  each  box  has  a  “this  side  up”  restriction  (e.g.,  Bischoff  and  Dows¬ 
land,  1982).  Even  without  the  “this  side  up”  restriction,  operational  considerations  may  dictate  the  use  of  ver¬ 
tical  layers  with  the  same  height.  Issues  of  stability  and  safety  of  the  boxes  imply  the  use  of  orthogonal  packing 
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(e.g.,  Dowsland,  1987a;  Nelissen,  1995;  Young-Gun  and  Maing-Kyu,  2001).  PLP  also  arises  in  some  cutting 
stock  and  floor  design  settings. 

Although  PLP  has  been  widely  studied  and  it  is  known  that  instances  can  be  divided  into  classes  with  the 
same  optimal  placement  pattern  (Dowsland,  1984),  no  procedure  to  generate  all  distinct  classes  for  a  given 
number  of  boxes  has  previously  been  reported.  Dowsland  (1987b)  works  with  a  subset  of  approximately 
8000  equivalence  classes  and  Scheithauer  and  Terno  (1996)  work  with  a  randomly-generated  subset  of  approx¬ 
imately  50,000  equivalence  classes.  A  common  set  of  restrictions  on  pallet  and  box  dimensions  first  proposed 
by  Dowsland  (1984)  has  been  used  by  other  authors  (e.g.,  Nelissen,  1993;  Scheithauer  and  Terno,  1996;  Mor- 
abito  and  Morales,  1998).  These  restrictions  are  on  the  aspect  ratio  of  the  pallet  (1  ^  X/Y^  2),  of  the  box 
( 1  <  a/b  <  4),  and  the  area  ratio  (1  <  (X*  Y)/(a  *  b)  <  51).  Nelissen  ( 1995)  and  Naujoks  (as  reported  by  Nelis¬ 
sen  (1995))  also  investigate  instances  where  51  <  (X  *  Y)l(a  *  b)  <  101  as  do  more  recent  papers  (e.g.,  Alvarez- 
Valdes  et  al.,  2005;  Birgin  et  al.,  2005;  Lins  et  ah,  2003)  that  apply  exact  algorithms  and  heuristics  to  a  com¬ 
mon  set  of  about  50,000  instances.  Recent  work  also  includes  detailed  analysis  of  upper  bounds  for  PLP 
(Letchford  and  Amaral,  2001). 

This  paper  defines  the  Minimum  Size  Instance  (MSI)  of  an  equivalence  class  of  PLP,  and  shows  that  every 
class  has  one  and  only  one  MSI.  We  develop  bounds  on  the  dimensions  of  box  and  pallet  in  the  MSI  of  each 
class.  Applying  our  newly-developed  bounds  on  the  MSI  dimensions,  we  present  an  algorithm  for  MSI  gen¬ 
eration  and  use  it  to  enumerate  all  3,080,730  equivalence  classes  with  an  area  ratio  (pallet  area  divided  by  box 
area)  smaller  than  101  boxes.  Previous  work  only  provides  bounds  on  the  ratio  of  box  dimensions  and  only 
considers  a  subset  of  all  classes:  this  limits  results.  Martins  (2003)  finds  all  instances  from  3,073,724  of  these 
3,080,730  classes  can  be  solved  easily.  Given  the  small  number  of  difficult  instances,  it  is  not  surprising  that 
many  have  been  previously  overlooked. 

2.  Efficient  partitions  and  equivalence  classes 

Some  PLP  instances,  with  different  dimensions,  possess  the  same  arrangement  of  boxes  in  an  optimal  solu¬ 
tion.  For  example,  the  arrangement  depicted  in  Fig.  1  is  an  optimal  solution  to  the  instance  (22,  16,  5,  3)  where 
the  shaded  regions  indicate  unused  (wasted)  areas  of  the  pallet.  The  same  arrangement  is  also  optimal,  for 
example,  to  instances  (30,  22,  7,  4)  and  (50,  36,  11,  7). 

Let  (n,m)  denote  an  ordered  pair  of  non-negative  integers  satisfying 

n*a  +  m*b^S  (1) 

for  a  pallet  dimension  S,  which  could  be  X  or  Y.  Such  an  ordered  pair  (n,m)  is  called  a  feasible  partition  of  S.  If 
n  and  m  also  satisfy 

0^S  —  n*a  —  m*b<b  (2) 


Fig.  1.  Optimal  solution  arrangement  for  instances  (22,  16,  5,  3),  (30,  22,  7,  4),  (50,  36,  11,  7),  and  all  instances  within  the  same  equivalence 
class. 
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then  (n,m)  is  called  an  efficient  partition  of  S  (Bischoff  and  Dowsland,  1982).  For  a  pallet  dimension  S,  the  set 
of  efficient  partitions  of  S,  denoted  by  E{S,a,b),  is  defined  to  be  the  set  of  all  feasible  partitions  (n,m)  satisfying 

n  e  {0, 1, . . . ,  [S/a|}  and  m  =  [(S  —  n  *  a)/b\. 

Dowsland  ( 1984)  shows  that  if  two  instances  of  PLP  possess  the  same  set  of  efficient  partitions  for  both  the 
pallet  width  and  length,  then  both  instances  share  the  same  set  of  optimal  solutions.  This  defines  a  relation  in 
the  set  of  instances  of  PLP,  which  is  reflexive,  symmetric,  and  transitive.  Therefore,  the  set  of  instances  of  PLP 
can  be  divided  into  equivalence  classes,  based  on  the  set  of  efficient  partitions.  If  a  solution  is  known  for  a  class 
representative,  then  this  solution  can  be  used  on  any  other  instance  in  the  class.  Because  multiplying  all  dimen¬ 
sions  by  an  integer  produces  a  new  instance  in  the  same  class,  it  is  easy  to  see  that  each  class  contains  infinitely 
many  instances. 

If,  in  addition,  n  and  m  satisfy 

n  *  a  +  m*  b  =  S, 

then  (n,m)  is  called  a  perfect  partition  of  S  (Dowsland,  1984).  In  general,  each  of  these  sets  can  be  empty,  but 
the  instance  of  an  equivalence  class  with  minimal  pallet  dimensions  contains  at  least  one  perfect  partition  for 
each  dimension,  X  and  Y  (Dowsland,  1984).  This  is  easily  observed  if  we  consider  an  arbitrary  instance  with¬ 
out  a  perfect  partition  for  a  given  dimension.  In  this  case,  we  can  reduce  the  corresponding  dimension  of  the 
pallet  without  altering  the  set  of  efficient  partitions.  This  implies  that  the  new  instance,  with  a  smaller  pallet 
dimension,  also  belongs  to  the  same  class. 

3.  Representing  equivalence  classes 

Because  instances  of  PLP  in  the  same  equivalence  class  share  the  same  set  of  optimal  solutions,  once  one 
instance  is  solved,  the  solution  can  be  stored  in  a  database  and  retrieved  whenever  a  solution  to  an  instance  of 
PLP  of  the  same  class  is  necessary  (Dowsland,  1987a).  Many  PLP  instances  are  easily  solved  so  storage  may 
only  be  necessary  for  difficult  instances.  The  most  straightforward  way  to  identify  an  equivalence  class  in  a 
database  is  to  encode  the  set  of  efficient  partitions  defining  the  class.  This  way,  given  a  new  instance,  it  is  pos¬ 
sible  to  compute  the  set  of  efficient  partitions  and  compare  it  with  the  entries  in  the  database.  One  possible 
problem  is  that  the  cardinality  of  this  set  increases  with  the  number  of  boxes  packed. 

Another  approach  is  to  select  a  unique  class  representative.  This  way,  only  four  integers  are  necessary  to 
represent  the  class,  independent  of  the  number  of  boxes  in  the  optimum  packing.  One  option  for  defining 
an  equivalence-class  representative  is  the  instance  that  minimizes  the  area  ratio,  the  Minimum  Area  Ratio 
Instance  {MARI).  But  the  minimization  problem  can  have  a  solution  at  an  open  boundary,  or  at  a  non-integral 
interior  point  (Dowsland,  1984).  In  these  cases  the  dimensions  of  the  MARI  can  only  be  approximated,  when 
using  integers.  Different  approximations  can  generate  different  instances  within  the  same  class,  complicating 
the  identification  process. 

Another  candidate  for  equivalence-class  representative  is  the  Minimum  Size  Instance  (MSI),  the  instance 
that  minimizes  the  dimensions  of  both  the  pallet  and  the  box.  We  say  (  X ,  Y ,  a,  b)  is  the  Minimum  Size  Instance 
of  a  class  if  for  all  instances  ( X ,  Y,a,b)  in  the  same  class,  X  <  X,  Y  ^  Y,  a  ^  a,  b  ^  b. 

4.  Existence  and  uniqueness  of  the  minimum  size  instance 

When  a  dimension  of  the  pallet  is  not  a  non-negative  integer  combination  of  the  box’s  dimensions,  Dows¬ 
land  (1984)  observes  that  the  dimension  of  the  pallet  can  be  reduced.  Let  G(S,a,b)  =  ma\{ij)^E(S^b) 
\i*a+j*b).  We  call  G(S,a,b )  the  Perfect  Partition  Equivalent  function.  Given  an  instance  ( X,Y,a,b )  of 
PLP,  the  reduced  dimensions  of  the  pallet  are  given  by  X*  =  G{X,a,b ),  and  Y*  =  G(Y,a,b).  Therefore,  if 
the  dimensions  of  the  box  in  the  MSI,  a  and  b,  are  known,  then  the  dimensions  of  the  pallet  are  given  by 
X  =  G(X,a,b )  and  Y  =  G(Y,a,b). 

We  show  that  the  MSI  is  unique  in  a  class  and  its  dimensions  can  be  easily  bounded,  simplifying  the  process 
of  enumerating  equivalence  classes. 
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Theorem  1.  Every  equivalence  class  of  PLP  has  one  and  only  one  MSI. 

Proof.  We  initially  show  that  there  is  no  more  than  one  MSI  in  each  class.  Then  we  show  that  every  class  has 
at  least  one  MSI. 

Suppose  (X\,Y\,a\,b\)  and  (X2,  Y2,a2,b2)  are  two  MSIs  in  an  equivalence  class.  By  definition,  both 
instances  minimize  all  dimensions  of  the  pallet  and  the  box  (X\  d-  X2,  Ih  Y2 ,  ct\  ^  $2,  ci\  a2  and  X2  st  X\, 
Y2  X  Y\,  a2  X  0\,  a2  <  ax),  implying  X,  =  X2,  Yx  =  Y2,  a\  —  a2,  rq  =  a2.  Therefore,  if  there  is  a  MSI,  it  is 
unique. 

Now  consider  an  equivalence  class.  Because  the  dimensions  of  the  pallet  in  the  MSI  are  a  function  of  the 
dimensions  of  the  box  in  the  MSI,  the  only  way  for  a  class  not  to  have  an  MSI  is  if  there  exists  one  instance, 
say  {X\,  Y\,a,b\),  with  minimum  length  for  the  box  (i.e.,  a  <  a  for  all  instances  (X,  Y,a,b)  in  the  same  class) 
and  another  instance,  (X2lY2,a2,b),  in  which  the  box  has  minimum  width  (i.e.,  b  b  for  all  instances 
(X,  Y,a,b)  in  the  class).  In  this  case,  a2>  a  and  b\  >  b.  The  strict  inequalities  hold  because  otherwise  at  least 
one  of  the  instances  would  have  the  box  with  both  minimum  dimensions.  As  both  instances  belong  to  the  same 
class,  E(X\,a,bi)  =  E(X2,a2,b)  and  E(Y\,a,b\)  =  E(Y2,a2,b).  We  show  that  the  MSI  can  be  identified  from 
these  two  instances. 

As  Dowsland  (1987a)  shows,  a  scaled  instance  of  PLP  remains  in  the  same  equivalence  class.  After 
scaling,  the  dimensions  of  the  pallet  and  box  may  no  longer  be  integers.  Normalizing  the  width  of  the  box 
to  1  in  the  above  instances,  we  obtain  instances  (X\/b\,Y\/b\,a/b\,\)  and  (X2/b,Y2/b,a2/b,\). 
Because  b\>b  and  b  >  0,  then  \/b  >  \/b\,  and  this  result  together  with  a2  >  a  give  us  a2/b  >  a/b  >  a/b\. 
Because  an  equivalence  class  is  a  convex  set  (Nelissen,  1993),  there  is  an  instance  (Xq,  Yu,  a/b.  1)  in  the  class. 
If  we  multiply  the  dimensions  by  b  we  obtain  the  instance  (X0  *b,Y0  *b,a,b).  We  can  apply  the 
perfect  partition  equivalent  function,  obtaining  X  =  G(X0  *  b, a , b)  =  maxp^^j,^,)!! *a+j*b}  because 
E(X o  *  b,  a,  b)  =  E{X\,a,  b\)  and  Y  =  G(To  *  b,  a,  b)  =  max(/,y)e£(r1,a,*i){i  *a+j*b}  because  E(Y o  *  b,  a,  b)  = 
E(Y\,a,  b\).  The  instance  ( X ,  Y,  a,  b)  satisfies  the  requirements  to  be  the  MSI  of  the  class.  Therefore,  the  class 
has  a  MSI.  □ 

5.  Bounds  on  the  dimensions  of  the  MSI  of  an  equivalence  class 

For  instance  (X,  Y,a,b),  let  Ax=  \X/a\,  Av=  |_ Y/a\ ,  Bx=  [X/b\,  and  By=  [ Y/b\ . 

Dowsland  (1987a)  shows  that  a  <  Bx  +  1  and  b  <  Ax  +  1  when  considering  the  set  of  ratios  a/b  correspond¬ 
ing  to  equivalence  classes.  While  these  limits  on  a  and  b  bound  the  ratio,  they  do  not  bound  a  and  b  in  an 
equivalence  class.  For  example,  instance  (104,90,  15,  13),  MSI  of  its  class,  where  Ax  =  6,  Aym  6,  Bx=  8, 
Br  =  6,  with  b  =  Ax  +  Ay+  1  and  a  —  Bx  +  By+  1.  Theorem  2  shows  these  are  upper  bounds  for  any  MSI. 

Theorem  2.  b  d  dx  +  Aj:  4  I  and  a  <  Bx  +  By  +  1 . 

Proof.  Given  instance  (X,  Y,a,b),  the  optimal  solution  to  the  integer  program  below  is  the  MSI  for  its  equiv¬ 
alence  class.  Using  the  optimal  solution  to  its  linear  programming  relaxation,  we  show  how  to  construct  a  PLP 
instance  from  the  equivalence  class  that  satisfies  the  bounds  of  Theorem  2. 

Indices 

i  efficient  partitions  on  length,  1  =  0,...,  Ax, 

f  efficient  partitions  on  width,  /=  0, . . . ,  Ay- 

Data 

pxj  number  of  boxes  with  their  largest  dimension  oriented,  vertically  in  partition  i  of  the  length, 

pyf  number  of  boxes  with  their  largest  dimension  oriented  horizontally  in  partition  /  of  the  width. 

Variables 

(X,  Y,a,b)  variables  for  (X,  Y,a,b). 
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Formulation 


Minimize  b 
subject  to 


X  —i*  a  —pXj  *  b  ^  0,  Vi  €  {0, 1, . . .  ,AX}, 
i  *  a  +  (pxi  +  1 )  *  b  —X  5=  1,  V/  e  {0,1,...  ,AX}, 
(Ax+l)*a-X  ^  1, 

Y-f*a  —pyf  *b^  0,  V/  G  {0, 1, . . .  ,Ay}, 
f  *a+(pyf+  1)  *  b-  Y  ^  1,  V/  e  {0,  l,...,Ay}, 
(Ay+l)*a-Y^  1. 

X,Y,  a,  b  integer 


(51) 

(52) 

(53) 

(54) 

(55) 

(56) 


We  know  from  Theorem  1  that  we  can  minimize  a  or  b  and  obtain  the  MSI,  here  we  minimize  b.  The  con¬ 
straint  sets  (SI)  and  (S4)  ensure  feasible  partitions,  Eq.  (1),  and  are  called  fitting  constraints.  The  other  con¬ 
straint  sets  ensure  efficient  partitions,  Eq.  (2),  and  are  called  efficiency  constraints.  We  call  the  linear 
programming  relaxation  (Primal). 

Inspection  of  constraints  (SI)  and  (S2)  reveals  that  the  addition  of  the  constraints  corresponding  to  the 
same  value  of  i  in  each  set  bounds  b  below  by  1,  i.e.,  ( X  —i  *  a  —pxt  *  b)  +  (i  *  a  +{pxi  +  1)  *  b  -X)  =  b  ^  1. 

Therefore,  when  (Primal)  is  feasible,  it  has  an  optimal  solution  (X*,  Y*,a*,b*),  with  objective  function  value  at 
least  1.  In  this  optimal  solution,  at  least  four  constraints  are  binding  because  it  is  a  four-dimensional  linear 
program  and  the  variables  have  no  non-negativity  constraints. 

Because  the  MSI  of  a  class  has  one  perfect  partition  in  each  dimension,  then,  at  least  two  fitting  constraints, 
one  in  the  length  (SI)  and  one  in  the  width  (S4),  are  binding  in  an  optimal  solution.  Also,  at  least  one  efficiency 
constraint  is  binding  -  from  (S2),  (S3),  (S5)  or  (S6). 

The  dual  of  (Primal)  has  only  four  rows  because  the  primal  has  only  four  variables.  This  makes  it  easier  to 
work  with  the  dual. 

Let  s  =  (sito,  ■  •  ■  ,51^,52.0,  ■  ■  •  ,S2,ax,S3,sA,o,  •  ■  ■  ,‘5'‘Mv,‘Wo,  •  -  -  ,J5^)J6)i  be  the  vector  of  dual  variables  of 
(Primal).  The  vector  shows  six  groups  of  variables  corresponding  to  the  first  and  last  dual  variable  for 
constraint  sets  (SI)  (S6).  The  dual  (Dual)  of  (Primal)  is  given  by 


Max 


Ay 


-  •?!/)  *  i  +  (Ax  +  1)  *  S3  +  ^2(s5f  ~  SAf)  *f  +  (Ay  +  1)  *  s6  =  0, 


Ay 


^(O,.  +  1)  *  s2i -pXi  *  Sii)  +  ^2{{pyf  +  1)  *  S5f  —  pyf  *  sAf)  =  1 


Ay 


subject  to 
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Let  P  be  the  matrix  of  technological  coefficients  of  (Primal). 

Then  PT 

has  the  following  structure: 

'  0 

■  -Ax  0 

Ax 

Ax  +  1 

0 

~Ay  0  . 

Ay  Ay  +  1 

-px  0  •  • 

•  ~PXAx  PX  0  +  1  ■ 

■  ■  PXAx  +  1 

0 

-pyQ  ■  ■  ■ 

-PVa,  P)’o  ■ 

■  •  PVa,  +  1  0 

1 

1  -1 

-1 

-1 

0 

0  0  . 

0  0 

0 

0  0 

0 

0 

1 

1  -1  . 

-1  -1 

In  every  optimal  dual  solution,  there  is  an  optimal  basis  that  contains  one  column  corresponding  to  a  perfect 
X-partition,  and  another  column  corresponding  to  a  perfect  7-partition.  This  follows  from  (Primal),  in  which 
there  is  always  one  binding  X-partition  row  and  one  binding  7-partition  row.  Therefore,  two  of  the  columns  in 
the  basis  look  like 


1  0 
.0  1. 

Also,  the  optimal  basis  contains  at  least  one  column  corresponding  to  a  binding  efficiency  constraint  in  the 
primal,  or  otherwise  the  dual  objective  function  would  have  value  zero.  Therefore,  the  basis  contains  at  least 
one  of  the  following  columns: 


-1 

or 

0 

.  0  . 

.-1. 

Considering  the  conditions  above,  the  optimal  basis  has  one  of  the  following  layouts,  not  considering  an  ex¬ 
change  in  the  last  two  rows: 


10-1-1 

.  (2) 

10-11 

,  (3) 

10-10 

,  or  (4) 

10-10 

.0  10  0  . 

1 

O 

o 

o 

1 _ 

1 

O 

O 

_ i 

1 

O 

o 

1 _ 

Let  B  be  a  basis,  and  let  sb  =  (sp,sg,sr,ss)  be  the  corresponding  vector  of  basic  variables.  The  coefficients  of 
the  first  two  rows  of  the  basis  are  identified  as  bkj,  k  =  1,2,  /  =  1,2, 3,4,  where  k  indicates  the  row  and  /  the 
column.  If  the  column  corresponds  to  a  fitting  constraint  in  (Primal),  its  coefficients  are  identified  by  bkj. 

The  cost  coefficients  for  the  (Dual)  objective  function  take  value  0  for  fitting  constraints,  and  1  for  efficiency 
constraints.  Therefore,  the  vectors  of  cost  coefficients  corresponding  to  the  cases  above  are  respectively 

(1)  (0,0, 1,1), 

(2)  (0,0, 1,0), 

(3)  (0,0, 1,0)  or 

(4)  (0,0, 1,1). 

In  all  four  cases,  we  can  reduce  the  system  of  equations  to  a  2  x  2  system. 

Case  (7).  The  system  is  given  by 

*1,1  *  Sp  +  *1,2  *  Sq  +  *!, 3  *  Sr  +  *1,4  *  ss  =  0, 

*2,1  *  Sp  +  *2,2  *  Sq  +  *2,3  *  Sr  +  *2,4  *  Ss  =  k 
Sp  —  Sr  —  =  0, 

Sq  =  0. 
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We  use  sq  —  0  and  sp  =  sr  +  ss,  and  obtain  the  following  system: 

(ft  13  +  fti  i)  *  Sr  +  (fti  4  +  fti  i)  *  Ss  =  0, 

(h2,3  +  ft2,l)  *  sr  +  (h2,4  +  h2, l)  *  Ss  =  1. 

Case  (2).  The  system  is  given  by 

ftl.l  *  $p  ~ t“  hi,2  *  Sq  T  hi, 3  *  Sr  ft  14  *  S,  =  0, 

ft2, 1  *  Sy,  +  ft2,2  *  S?  +  ft2)3  *  Sr  +  ft2i4  *  Ss  =  1 , 

Sp  Sy  ~\~  S s  0, 

Sq  =  0. 

We  use  Sq  —  0  and  sr  =  sp  +  ss,  and  obtain  the  following  system: 

(hi, 3  +  ftu)  *  sp  +  (hi, 3  +  hi, 4)  *  ss  =  0, 

(h2, 3  +  h2,i)  *  sp  +  (ft2, 3  +  ft2,4)  *  ss  =  1. 

Case  (5).  The  system  is  given  by 

fti.i  *  sp  +  ft i,2  *  sq  +  hi, 3  *  sr  +  ft  1,4  *  ss  =  0, 

ft2, 1  *  Sy,  +  ft2,2  *  S?  +  ft2,3  *  Sr  +  ft2,4  *  Ss  =  1 , 

Sy,  Sy  0, 

Sq  T  Ss  —  0. 

We  use  sp  =  sr  and  sq  =  —ss,  and  obtain  the  following  system: 

(hi, 3  +  hi,i)  *  sr  +  (hi, 4  -  hi,2)  *  ss  =  0, 

(h2, 3  +  h2,i)  *  sr  +  (h2,4  —  h2,2)  *  ss  =  1. 

Case  (4).  The  system  is  given  by 

hi  1  *  Sp  +  fti,2  *  s^  -T  hi, 3  *  sr  -T  hi  4  *  ss  =  0, 

h2,i  *  sp  +  h2,2  *  s?  +  h2,3  *  s,.  +  h2, 4  *  ss  =  1 , 

Sy,  S;-  -  0, 

s^  —  0. 

We  use  sp  —  sr  and  sq  —  ss,  and  obtain  the  following  system: 

(hi, 3  +  hi,i)  *  sr  +  (hi, 4  +  hi,2)  *  ss  =  0, 

(h2, 3  +  h2, 1)  *  sr  +  (h2, 4  +  ft2,2)  *  ss  =  1. 

In  each  of  the  four  cases,  let  d  be  the  determinant  of  the  2x2  matrix.  The  matrix  is  a  basis  and  all  elements  are 
integers,  so  \d\  >  1.  Therefore,  the  respective  optimal  objective  function  values  are 

(1)  b'  =  sr  +  ss  =  ( — (hi.i  +  hi, 4)  +  (hi,i  +  hi,3))/d  ^  |(hi,3  —  hi.4)/ c/|  ^  (Ax  +  l)/|<f|, 

(2)  b*=sr  =  (hj  1  -  hi,4)/d  ^  |(hj,i  -  fti,4)/rf|  <  Ax/\d\, 

(3)  b*  =  sr  =  (hi,2  -  fti,4)/d  ^  |(hi,2  -  hi,4)/rf|  ^  Ax/\d\,  and 

(4)  ft*  =  sr  +  ss  =  ( — (hi .4  +  hi,2)  +  (hi, 3  +  b\  \))/d  ^  |((hi,3  +  hi,i)  —  (hi, 4  +  hi,2))/d|  ^  (Ax  +Ay  +  l)/|rf|. 

Because  all  elements  of  the  basis  are  integers,  h*  |d|  is  an  integer. 

In  case  (1),  both  coefficients  are  non-negative,  corresponding  to  efficiency  constraints,  and  one  can  take 
value  0.  Therefore,  the  maximum  value  of  b*  \d\  is  Ax  +  1. 

In  cases  (2)  and  (3),  all  coefficients  are  non-positive,  corresponding  to  fitting  constraints,  and  one 
coefficient,  in  each  case,  can  take  value  0.  Therefore,  the  maximum  value  of  b*  \d\  is  Ax. 
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In  case  (4),  if  we  consider  the  order  of  the  rows  as  presented,  the  result  follows  because  (£1,3  + hi, 1)  is 
computed  over  constraints  in  X,  taking  maximum  value  Ax  +  1,  and  (h] .4  +  hi, 2)  is  computed  over  constraints 
in  Y,  taking  minimum  value  —Av.  If  the  order  of  the  rows  is  exchanged,  the  values  obtained  are  also 
exchanged,  with  (hi, 3  +  b  1,1)  taking  maximum  value  Av  +  1  and  (hi, 4  +  hi,2)  taking  minimum  value  —Ax. 
Considering  all  four  cases,  we  verify  that  Ax+  Ay  +  1  is  an  upper  bound  on  b*  \d\.  Scaling  the  optimal 

(Primal)  solution  by  \d\,  \d\{X* ,  Y*,a*,b*)  is  anjnteger  solution  that  is  feasible  to  (Primal).  We  therefore  have 
an  instance  in  the  equivalence  class  with  b  =  b*  \d\  and,  by  definition,  h  <  h  so 

h  ^  Ax  +  Ay  +  1 .  (3) 

If  we  minimize  a  in  (Primal),  then  the  upper  bound  for  a  is 

a  ^  Bx  +  Bv  +  1 .  D  (4) 

6.  Identifying  the  MSI 

Given  an  instance  of  PLP,  {X,  Y,  a,  b)  or  the  corresponding  set  of  efficient  partitions,  we  test  values  for  h, 
starting  at  1,  and  compute  the  other  variables,  until  the  MSI  is  found.  Our  algorithm  operates  with  two  main 
loops.  The  outer  loop  selects  values  for  h,  from  1  to  min  jh.  Ax  +  Av  +  1 ) .  For  each  value  of  h,  we  compute  the 
range  of  a,  R-a,  given  by  R-a  =  {a  €  Z+  :  |  Bx  *  b/(Ax  +  1)J  <  a  ^  \[{BX  +  1)  *  &  —  1  )/Ax~\}.  The  second  loop  enu¬ 
merates  values  for  a  in  R„.  With  a  and  h,  we  compute  X  =  G(X .  a,  b).  If  the  selected  values  for  X,  a,  and  h 
satisfy  the  efficiency  constraints^  we  repeat  the  same  computations  for  Y,  or  otherwise  try  the  next  value 
for  a.  If  the  inequalities  for  Y  are  not  satisfied,  we  continue  with  the  procedure:  otherwise,  instance 
{X ,  Y,a,b)  is  the  MSI.  This  algorithm  is  simple  to  implement,  with  time  complexity  Of/!4). 

7.  Generating  equivalence  classes 

We  enumerate  the  MSI  of  all  equivalence  classes  with  an  area  ratio  smaller  than  101  boxes  per  pallet.  Also 
we  use  the  MSI  to  uniquely  identify  each  class  and,  therefore,  record  only  one  instance  per  class. 

If  N  is  the  maximum  number  of  boxes  that  can  be  packed  on  a  pallet,  we  have 


A  x  +  Ay  +  1 

(5) 

Bx  +  By  ^  2  N. 

(6) 

We  recall  some  definitions  to  demonstrate  these  bounds.  Any  optimal  must  have  at  least  Ax  (Ay)  boxes  placed 
side  by  side  across  the  length  (width)  of  the  pallet,  so  Ax  *  Ay  <  N.  If  Ay  =  0,  then  Ax  <  N  and  Ax  +  Ay  <  N.  If 
Ay  ^  1,  then  Ax  C  N/Av  and  Ax  +  A  v  N/Av  +  Ay  <  N+  1.  Also,  Bx  ( By)  is  the  maximum  number  of  boxes 
that  can  be  placed  side  by  side  across  the  length  (width)  of  the  pallet.  Therefore,  Bv  A  Bx  N  and 

By  +  Bx  <  2  N. 

Our  PLP  Equivalence  Class  Generation  Algorithm  (PLP-ECGA)  has  six  main  loops,  and  uses  a  list  ordered 
lexicographically  by  b,  a  and  Y  to  maintain  the  distinct  generated  classes  for  given  values  of  a  and  b.  The  out¬ 
ermost  loop  determines  the  values  for  b,  from  1  to  N  +  2,  because  b  ^  Ax  +  Ay  +  1  (3)  and  Ax  +  Ay  <  N  +  1  (5). 
The  second  loop  selects  values  for  u,  from  b  +  1  to  2 N  +  1,  because  a  <  Bx  +  By  +  1  (4)  and  Bx  +  By  <  2 N  (6), 
except  when  b  equals  1 ,  when  a  can  also  be  equal  to  1 .  If  the  greatest  common  divisor  of  a  and  b  is  greater  than  1 , 
then  the  second  loop  proceeds  to  the  next  value  for  a.  Otherwise,  the  ordered  list,  with  all  generated  classes,  is 
emptied.  The  third  and  fourth  loops  select  among  all  the  possible  perfect  partitions  of  the  width  for  candidate  Y. 
The  fifth  and  sixth  loops  select  among  the  perfect  partitions  of  the  length  for  candidate  X.  If  instance  ( X ,  Y,  a,  b) 
has  an  area  ratio  bound  not  exceeding  N  and  has  not  been  generated  before  then  it  is  recorded  in  the  list. 

It  is  possible  to  verify  through  the  algorithm  that  the  number  of  equivalence  classes  is  bounded  by  a  poly¬ 
nomial  in  N,  albeit  a  large  polynomial.  There  are  0(A2)  ways  of  assigning  values  to  a  and  b,  corresponding  to 
the  number  of  pairs  of  relatively  prime  numbers  less  than  or  equal  to  2 N.  More  precisely,  the  number  is  given 
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Table  1 

Number  of  equivalence  classes  in  groups  based  on  the  maximum  number  of  boxes 


Maximum  number  of  boxes  in  each  group  Number  of  classes  in  each  group  Run  time  to  generate  all  classes  (seconds) 


10 

662 

0 

20 

7309 

3 

50 

216,095 

295 

100 

3,080,730 

14,667 

Table  2 

Distribution  of  values  of  h  in  the  MSI  in  each  class 


Number  of  boxes 

Number  of  classes 

6  =  1 

b^2 

b^5 

b^  10 

b^  20 

b^  50 

10 

662 

92 

276 

609 

662 

662 

662 

20 

7309 

520 

1760 

4873 

6659 

7309 

7309 

50 

216,095 

6362 

23,270 

71,686 

119,298 

182,870 

216,095 

100 

3,080,730 

46,300 

174,177 

544,004 

964,673 

1,710,574 

2,822,767 

by  Y^k= where  <j>{k)  is  the  Euler  phi-function,  which  gives  the  number  of  integers  less  than  k  that  are 
relatively  prime  to  k  (Gallian,  1998).  The  loops  corresponding  to  the  width  are  executed  0(  A-2)  times  for  each 
pair  a  and  b.  The  same  happens  with  the  loops  corresponding  to  the  length.  Therefore,  the  number  of  equiv¬ 
alence  classes,  with  area  ratio  bound  up  to  N,  is  bounded  above  by  a  sixth-degree  polynomial  in  N. 

The  instances  generated  with  the  PLP-ECGA  procedure  are  divided  in  groups  of  up  to  10,  20,  50  and  100 
boxes  per  pallet,  as  defined  by  the  area  ratio  bound.  Table  1  presents  in  the  second  column  the  number  of 
equivalence  classes  of  PLP  in  each  group.  The  third  column  contains  the  time  required,  in  seconds,  to  generate 
the  equivalence  classes  within  each  group  on  a  Pentium  III  600  MHz  personal  computer. 

Table  2  presents  the  distribution  of  classes  in  each  group,  where  the  MSI  is  defined  with  b  smaller  than  or 
equal  to  1,2,  5,  10,  20,  and  50.  The  first  column  defines  the  maximum  number  of  boxes  that  can  be  packed  in 
an  instance  in  the  group  of  classes  covered,  as  given  by  the  area  ratio  bound.  The  second  column  lists  the  num¬ 
ber  of  distinct  classes  in  each  group.  The  following  columns  present  the  number  of  classes. 

For  PLP  instances  with  area  ratio  smaller  than  101  boxes,  91%  of  the  equivalence  classes  have  a  MSI  where 
the  value  of  b  is  less  than  or  equal  to  50.  The  complete  set  of  instances  can  be  accessed  at  http://www.pallet- 
loading.org.  More  details  can  be  found  in  (Martins,  2003). 

8.  Conclusions 

In  this  paper,  we  define  the  Minimum  Size  Instance  (MSI)  of  an  equivalence  class  of  PLP,  and  show  that 
every  class  has  one  and  only  one  MSI.  This  makes  the  MSI  helpful  in  distinguishing  equivalence  classes.  We 
also  develop  bounds  on  the  dimensions  of  the  box  and  pallet  in  the  MSI  of  a  class.  Previous  work  only  pro¬ 
vides  bounds  on  the  ratio  of  box  dimensions.  Applying  the  newly  developed  bounds  to  the  MSI,  we  enumerate 
the  MSI  of  all  equivalence  classes  with  area  ratio  smaller  than  101  and  provide  some  statistics. 

References 


Alvarez-Valdes,  R.,  Parreno,  F.,  Tamarit,  J.M.,  2005.  A  branch-and-cut  algorithm  for  the  pallet  loading  problem.  Computers  & 
Operations  Research  32,  3007-3029. 

Birgin,  E.G.,  Morabito,  R.,  Nishihara,  F.H.,  2005.  A  note  on  an  L-approach  for  solving  the  manufacturer’s  pallet  loading  problem. 
Journal  of  the  Operational  Research  Society  56,  1448-1451. 

Bischoff,  E.E.,  Dowsland,  W.B.,  1982.  An  application  of  the  micro  to  product  design  and  distribution.  Journal  of  the  Operational 
Research  Society  33,  271-281. 

Dowsland,  K.A.,  1984.  The  three-dimensional  pallet  chart:  An  analysis  of  the  factors  affecting  the  set  of  feasible  layouts  for  a  class  of  two- 
dimensional  packing  problems.  Journal  of  the  Operational  Research  Society  35,  895-905. 

Dowsland,  K.A.,  1987a.  A  combined  data-base  and  algorithmic  approach  to  the  pallet-loading.  Journal  of  the  Operations  Research 
Society  38,  341-345. 


26 


G.H.A.  Martins,  R.F.  Dell  /  European  Journal  of  Operational  Research  179  (2007)  17—26 


Dowsland,  K.A.,  1987b.  An  exact  algorithm  for  the  pallet  loading  problem.  European  Journal  of  Operational  Research  31,  78-84. 

Gallian,  J.A..  1998.  Contemporary  Abstract  Algebra.  Houghton  Mifflin  Company,  Boston. 

Letchford,  A.N.,  Amaral,  A.,  2001.  Analysis  of  upper  bounds  for  the  pallet  loading  problem.  European  Journal  of  Operational  Research 
132,  582-593. 

Lins,  L.,  Lins,  S.,  Morabito,  R.,  2003.  An  L-approach  for  Packing  (/,  w)-rectangles  into  rectangular  and  L-shaped  pieces.  Journal  of  the 
Operational  Research  Society  54,  777-789. 

Martins,  G.,  2003.  Packing  in  two  and  three  dimensions.  Ph.D.  Dissertation,  Naval  Postgraduate  School,  Monterey,  California,  June. 

Morabito,  R.,  Morales,  S.,  1998.  A  simple  and  effective  recursive  procedure  for  the  manufacturer’s  pallet  loading  problem.  Journal  of  the 
Operational  Research  Society  49,  819-828. 

Nelissen,  J.,  1993.  New  approaches  to  the  pallet  loading  problem.  Working  paper,  RWTH  Aachen.  Available  from:  <ftp:// 
ftp.informatik. rwth-aachen.de/pub/reports/pallet. ps.Z>. 

Nelissen,  J.,  1995.  How  to  use  structural  constraints  to  compute  an  upper  bound  for  the  pallet  loading  problem.  European  Journal  of 
Operational  Research  84,  662-680. 

Scheithauer,  G.,  Terno,  J.,  1996.  The  G4-heuristic  for  the  pallet  loading  problem.  Journal  of  the  Operational  Research  Society  47,  51 1 — 
522. 

Young-Gun,  G.,  Maing-Kyu,  K.,  2001.  A  fast  algorithm  for  two-dimensional  pallet  loading  problems  of  large  size.  European  Journal  of 
Operational  Research  134,  193-202. 


